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Using a simple diagram technique we derive the electric-dipole echo amplitude from two-level 
systems with a quadrupole nuclear moment in glasses in an external magnetic field. We show, that 
due to the quadrupole moment interaction of a tunneling particle with a gradient of an internal 
electric field, the echo amplitude experiences oscillations in rather weak magnetic fields. With an 
increase of the magnetic field, when the Zeeman energy becomes larger than the quadrupole energy 
splitting, the average echo amplitude increases and saturates (for high magnetic fields) at some level 
which is above the average level of echo oscillations for small magnetic fields. 
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I. INTRODUCTION 

Recently it was experimentally observed that the electric-dipole echo amplitude in non-magnetic glasses oscillates 
as a function of applied magnetic field-' 2 . Similar behavior was found in insulating crystals with tunneling impurities 3 . 
In a recent paper 4 such behavior was attributed to quadrupole nuclear moments of tunneling particles interacting 
with the magnetic field and with the gradient of an internal electric field. The purpose of this paper is to formulate 
a general microscopic theory of this phenomenon taking into account a multi-level structure of tunneling systems 
in glasses coupled with a quadrupole nuclear moment. To this end we have developed a simple diagram technique 
describing the coherent echo signal. Making use of this technique, one can easily generalize known results for the echo 
amplitude from the usual two-level systems (TLS) to more general multi-level systems. 

II. DIAGRAM REPRESENTATION OF THE ECHO 

In this section we give a simple diagram representation of the two-pulse echo signal from an arbitrary multi-level 
system. We then apply this technique to a multi- level system formed by a TLS with a quadrupole nuclear moment. 
The microscopic details of the TLS-quadrupole interaction in glasses are discussed further down in the next section. 
First we start with a simple example of an echo in an ensemble of TLS. To avoid unnecessary complications we will 
use a simple perturbative approach regarding the applied electric field, similar to the one considered in Ref£. However 
all results can easily be generalized to strong electric fields. 

A. Echo in an Ensemble of TLS 

The wave function of a TLS in an external ac-electric field is a linear combination of the wave functions ip% and if 2 
for each level 

* = C m + C 2¥ > 2 , \d\ 2 + \C 2 \ 2 = 1. (1) 

Prior to the action of the first electric pulse we have C\ = 1 and C 2 = 0. Then in the electric field the time variation 
^ ■ of the amplitudes C\ and C 2 obeys the equations 
5_i ' 

ih^ = E 1 C 1 + V(t)C 2 , ih^-=E 2 C 2 + V(t)C 1 . (2) 
at at 

For the off-diagonal transition matrix element V(t) we have a following expression during the electric pulse 

V(t) = V, 2 coswt, where V h2 = (F lj2 • m)^. (3) 

Here Fi j2 is the electric field amplitude of the first or the second electric pulse, respectively, m is the dipole moment 
of the TLS, A the tunneling splitting, E = E 2 — Ei the TLS energy, and E\ and E 2 are the energies of the ground 
and the excited states of the TLS, respectively. The electric field frequency is assumed as fiuj w E. To simplify our 
equations we will put h = 1. 
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After the action of the 1st electric pulse C 2 acquires in the first approximation a finite value proportional to 
amplitude V\ which is assumed to be small. During the time interval between the first and second electric pulses 
(0 < t < T21) we have 

C 2 oc {-iVx) -e' tE '\ & « 1 • e~ tElt - 0{V?). (4) 
In Eq. J2J for C 2 we have omitted the resonance factor 

^2i(n) = - sin ^ exp ^-iE 2Tl + i^-j , (5) 

where r% is a duration of the first pulse, and z — E — lj is a detuning from the resonance. For a transition between the 
second and the first level (under the action of the second pulse, see below) the corresponding resonance factor reads 

/3i2(t 2 ) = -- sin ^ exp (-iEtT 2 - ■ (6) 

These resonance factors are very important in describing the form of the echo envelope^. But in the present paper 
we are mainly interested in the echo amplitude as a function of the time delay 721 between the two pulses. To make 
our equations as simple as possible, we assume that t 2 \ is much larger than the duration of any pumping pulse, r,-. As 
we will show below, if the electric pulses Tj are sufficiently short (the exact criterion for a TLS-quadrupole multi-level 
system will be given in the last section) such resonant factors are not important for our purposes and we will not take 
them into account in the following analysis (though they can easily be included if necessary). 

Let us now discuss what happens just after the action of the second electric pulse V 2 cosu>t. Using Eq. (J2J we 
calculate the corresponding variations of C\ and C 2 ■ 

8d cx (-iV 2 ) • C 2 (r 21 ) cx (-iVa) ■ (-iV 1 )e~ iE ^, 

(7) 

SC 2 cx {-iV 2 )-C 1 {T 21 )cx{-iV 2 )-e- iE ^\ 

(here for the reasons mentioned above we have omitted the resonance factors /3i 2 (r 2 ) and /3 2 i(t 2 ) in 8C\ and SC 2 , 
respectively). Finally, at some moment t after the 2nd pulse we have 

SC-l cx -V{V 2 e- iE ™ ■ e -^i(*-^i) ; 5 c 2 « -iV 2 e~ lElT21 ■ e -<^(*-^0. ( 8 ) 

The amplitude of the two-pulse echo from one TLS is determined by the average value of the product 8C\8C 2 {i.e. 
by the off-diagonal density matrix element) which is proportional to 

p ccho cx 5d5C* 2 cx -iVxVie^-^™ ■ e i(s 2 -B 1 )(t-r 21 )_ ^ 

Summing over all resonant TLS we have for the two-pulse echo amplitude 

Pcchocx-z^^e^- 2 ^). (10) 

TLS 

For t = 2r 2 i the amplitude of the echo has a sharp maximum, P e cho cx —iV\V 2 N, where N is the number of resonant 
TLS (i.e. TLS with fruj sa E). 

Let us now give a simple diagram representation of this result. It will allow us to generalize our approach to the more 
complicated case of an echo in a multi-level system formed, for example, by a tunneling particle with a quadrupole 
nuclear moment. For simplicity, we assume that the electric field amplitudes F 12 of the two electric pulses are parallel 
to each other, i.e. Fj = i^e, where e is a unit polarization vector. Then we have the following expression for the 
induced dipole moment echo amplitude from a single TLS 

Pecho cx {-iF x )a 21 ■ e - iE ™ ■ (-iF 2 )a 12 ■ e -^(*-«0 . \ e -iBim . (-iF 2 )a 21 ■ e -^(*-«0] * . a2l = 

(11) 

= -iF^a^a^e 1 ^-^-^-^-^^ = -iF^a^a^e 1 ^- 2 ^ . 

Here otij are the off-diagonal dipole transition matrix elements projections on the direction of the electric field. In the 
considered case they are real quantities and given by the usual expression 

«12 = a 2 i = (e- m )-gT- (12) 
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FIG. 1: A diagram for a two-pulse echo from a TLS. 




FIG. 2: A diagram for the two-pulse echo from TLS-quadrupole multi-level system. 



Generally they obey the relations a\2 = a 2 i- 

Eq. Ijlljl is similar to Eq. and can be depicted by the diagram shown in Fig. ^ Each full or dashed vertical 
line in the diagram corresponds to a transition from level i to level k (i,k = 1,2) under the action of the first or 
second electric pulse. To each vertical line we ascribe a factor (—iFi, 2 )aki (depending of the pulse). For example, 
the line DA corresponds to the factor (— iF\)ct2i and the full line BE to the factor (— zi^o^- The dashed line EB 
corresponds to the factor (— iF 2 )a 2 \ ■ Finally the dotted line FC (the echo signal) corresponds to the factor a 2 i- 

Each horizontal line in this diagram (full or dashed) corresponds to free TLS dynamics in the time intervals 
between or after the pulses. Therefore, an appropriate exponential factor is ascribed to each such line. The factor 
exp(— iE2T2i) corresponds to the line AB and the factor exp(— iE\T 2 \) to line DE. The factors exp[— iE 2 (t — r 2 x)] 
and expj— iE\{t — 121)] are represented by the lines BC and EF, respectively. In the final expression all factors 
corresponding to dashed lines should be taken as complex conjugated. It is easy to see that all full lines (vertical and 
horizontal) contribute to the amplitude 8C\ and all dashed lines to SC2- 



B. Echo in a Multi-Level System 



Using these simple diagram rules we can now easily find the contributions to the two-pulse echo signal from a 
multi-level system. Let us, for example, consider a multi-level system as shown in Fig. [21 It consists of two identical 
groups of N levels shifted vertically against each other by some energy E (playing the role of the usual TLS energy) . 
Inside the groups the positions of the levels are arbitrary, i.e. they are not necessarily equidistant. As we will see in 
the next section, such a multi- level system describes a TLS with a quadrupolc nuclear moment. 

According to the rules formulated above, the diagram in Fig. [21 gives the following partial contributions to the 
two-pulse echo signal 



5rii,fcij2,i2(t) = -4^{-iFi)a j2 ,a ■ e^ E ^ ■ {-iF 2 )a kl>j2 ■ e -«^(*-^) x 



1 p~i-E.lT,! 



N 



(-iF 2 )a l2ti i ■ e- iE ^- T ^ 



«12,kl 



(13) 



= -jjFiF%a j2 ,n ■ a k i,j2 ■ a* l2 il ■ a/ 2 ,fci • exp [i(Ei 2 - E kl )(t - r 2 i) - i(E j2 - Enfoi] 
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The factors 1/y/N in this formula correspond to the case when the low-energy group levels are equally populated 
(with a probability l/N) and the high-energy group levels are empty. This is the case for low enough temperatures, 
T <C E. On the other hand, to satisfy the previous conditions the temperature must be much larger than the width 
of the energy splitting in the groups (this width is of the order of the quadrupole splitting AEq), i.e. T ^> AEq. The 
latter inequality corresponds to the usual experimental situation. The two conditions are compatible if E ^> AEq. 
In the usual experiments this inequality is obeyed since the resonance frequency hto E 3> AEq. The limitation 
E T is not a crucial. The final result can be easily generalized to the case T ~ E by including thermal occupation 
numbers. 

Since the two group of levels are identical and only shifted by the TLS energy E we have E12 = E + En and 
Ej2 = E + Eji. Then Eq. l(T3|l can be rewritten as 

7r«, fc ij2,i2(t) = -^FiFZocfpa^a^af^ • exp [iE(t - 2r 21 ) + i(E lx - E kx ){t - r 21 ) - i(E n - E %1 )t 21 ] . (14) 

For convenience, we write the indices (1,2) of the level groups as superscripts. From Eq. 114(1 it follows, in the case 
when E 3> AEq, that the echo signal appears at t — 2t2\. At this time, summing over all possible combinations of 

(12) *(21) 

vertical transitions (i, j, k, I) between different levels and taking into account that a\j = , the total contribution 

to the echo signal from one multi-level system ism 



Pccho(2r 2 i) oc -— FiF 2 2 ^2 



e i(Ei-E k )T 21 



,k 



E 



(12) *(12) iE jT21 



(15) 



This expression differs from the similar one, Eq. (8) of Ref. 4 



III. TLS INTERACTION WITH NUCLEAR QUADRUPOLES 

In this section we consider a TLS with a nuclear quadrupole electric moment in external electric and magnetic 
fields. First we derive a Hamiltonian describing the interaction and then will apply the perturbation theory approach 
in respect to this interaction to describe the echo phenomenon in this system. 

A. General Relations 

The Hamiltonian of a nuclear electric quadrupole interacting with a gradient of an internal electric field has the 
usual form& 

= Qik ¥>ik, ¥>ik = Q r Q r S = Atp = 0. (16) 



Here tp is an electrostatic potential at the nuclear site and the traceless tensor 



Qik 



6.7(2.7-1) 



3 - - 

-{J l J k + JkJi) - So, J (J + 1) 



(17) 



is the operator of the nuclear quadrupole electric moment, J is the operator of the nuclear spin. 
The interaction of a nuclear magnetic moment with a magnetic field H can be written as 

Hn = -jhJ H = M H (18) 

where 7 is the nuclear gyromagnetic ratio and M is the nuclear magnetic moment operator. Both Hamiltonians lltjl) 
and l|18fl arc Hermitian N x N matrices, with N = 2 J + 1. 

The total Hamiltonian of a tunneling particle (or a group of particles) with a quadrupole nuclear moment in an 
external electric and magnetic field can be written in as 

Wtot = V(x) - F • d — + Qik tp ik (x) - M • H. (19) 

Xo 

Here x is a generalized coordinate of the tunneling particle, V(x) is a soft atomic double- well potential 9 , F is the 
applied electric field, dorr/a^ is the particle electric dipole moment, and xq ~ lA is of the order of interatomic 
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distance. The internal electric field gradient tensor at the site of the nucleus, (pik(x), is a function of the generalized 
particle coordinate x. 

Since the relative displacement of a tunneling particle x/xq <C 1 we can expand <p%k{x) in the Taylor series and 
limit ourselves to the linear approximation 

¥>i*(x) = Wfe(0)+¥/ tt (0)— ■ (20) 

Xq 

The second rank tensors <^zfe(0) and ip' ik (0) are independent of each other and are of the same order of the magnitude, 
ip' ik (0) ~ <fiik(0). Since the electric field gradient is a traceless tensor, <pu(x) = for any x, the same property holds 
for the tensors <fiik(0) and ip' ik (0), i.e. <pu(0) — <p'u(0) = 0. 
As a result the total Hamiltonian Ijl9(l becomes 

Htot = V(x) - M ■ H + Q lk ip ik (0) - F • d — + Q lk <p' ik (0) — . (21) 

Xo Xo 

The first term in this expression is the potential energy of the tunneling particle. The second and the third terms 
describe the interaction of the nuclear magnetic moment with the magnetic field and of the quadrupole nuclear 
moment with the average internal electric field gradient, respectively. The fourth term describes the interaction of the 
particle with the external electric field and finally the last, and for our theory most important term, accounts for the 
interaction of the quadrupole nuclear moment with the particle "orbital" motion in the soft atomic potential V(x). 
We will see that this last term is responsible for the electric-dipole echo oscillations in a magnetic field H. 

B. TLS Approximation 

To proceed further we will use the usual TLS approximation, keeping in mind that the particle moves in a nearly 
symmetric double-well potential with two minima at x\ and X2 ~ —x\. Taking the zero of the potential energy V(x) 
as the average of the values V(x\) and Vfez), we can write 

V(x 1 ) = ±A, V(x 2 ) = -±A, (22) 

where A is the energy difference between the two minima (the TLS asymmetry). Taking the tunneling under the 
barrier into account, we can substitute in this approximation the double-well potential energy V(x) by a 2 x 2 matrix 

where Ao is the usual tunneling amplitude. 
Similarly way we can write 



k,|( I _ j ] kik,- 124) 



where a, is the Pauli matrix. As a result 



-F-d — =^F-m£ z , m = d ^, (25) 
Xo x 

where m is the TLS dipole moment and 

Q ik ^ fc (0) — =► a z (8 V Q , Vq = -Q lk ^.(0) l^ii. (26) 
xo xo 

Introducing the notation 

Wq = -M • H + Qik ^ife(O), (27) 
we get for the Hamiltonian H21J1 in the TLS approximation 

Utot = -( J ®lQ + l a ®W Q + F ■ md z ®l Q + d z ® V Q . (28) 
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Here 1q is a N x N unit matrix in the space of nuclear spin J and 1 CT is a 2 x 2 unit matrix in the TLS, a space. The 

symbol (g) indicates the direct product of two matrices. As a result 7itot is a 27V x 2N Hermitian matrix in a conjoint 
space. The last term in this equation describes the interaction of the nuclear quadrupole with the TLS motion. It 
will be responsible for the echo oscillations in a magnetic field. In the following analysis we will consider this term to 
be small and treat it by standard perturbation theory. 

To proceed further, we should diagonalize the first two terms in Eq. I|28|l . The first term can be diagonalized using 
a standard unitary transformation in TLS space 



A 

-A 



-A 
-A 



s-' = 



E 
-E 



S„ = 



cos 9 — sin 9 
sin 9 cos 9 



(29) 



where E = ^Ag + A 2 and tan 20 = A /A, sin 29 = Aq/E. Under this transformation the a z matrix in the third and 
the fourth terms in Eq. (|28[) is transformed, as usual, to 



S(j<JzS„ — S a ( q . ) s„ 



( A Ao \ 

E E 

A A 
~Ej 



(30) 



In a similar way we can diagonalize the second term in Eq. (|28|l and transform the fourth term as follows 



SqWqSq^Wq, and S Q V Q S? 



Vq- 



(31) 



Here Wq is a diagonal matrix in the nuclear spin space N x N. It gives the nuclear quadrupole energies in the 
average (over the two minima) internal electric field gradient and in the external magnetic field H. The unitary 
transformation Sq can be found in general only by numerical diagonalization of the matrix Wq. The transformed 

matrix Vq describes the interaction of the TLS with the nuclear quadrupole moment. 

Finally as a result of these two independent unitary transformations the total TLS-quadrupolc Hamiltonian reads 



^tot — ~ 



E 
—E 







f * 


A \ 




t * 


A \ 








E 


E 




E 


E 






+ T CT <g> W Q + F • m 


A 


A 




A 


A 


®Vq 








~e) 






~e) 





f 



V 



(32) 



Here the first two terms are diagonal matrices. Together, they give two identical groups of levels (determined by 

the eigenvalues of Wq) shifted from one another by the TLS energy E. For convenience we have introduced the 
abbrevations / and V for the last two terms. This notation will be used in the next section. 



IV. DIPOLE TRANSITION MATRIX ELEMENTS 



According to Eq. (|15|) the electric-dipole echo amplitude is determined by the off-diagonal dipole transition matrix 

(12) 

elements, a ik , between the levels in the TLS-quadrupole multi- level system. In the present section, we will calculate, 
using the Hamiltonian l|32() . these matrix elements by standard perturbation theory. The perturbation will be the 
last term in Eq. V. 

Let us consider the TLS-quadrupole multi-level system shown in Fig. It consists of two identical groups of levels 

i 12) 

(1) and (2), shifted by the TLS energy E. We are interested in the off-diagonal transition matrix elements /Am 
between these two groups of levels. In zero order (V = 0) transitions are induced by the third term in Eq. I|32|) . 
However, the only non-zero off-diagonal matrix elements are the ones for transitions between identical levels in the 
two groups 

flT 2) = F • m^. (33) 
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FIG. 3: A TLS-quadrupole multi-level system. 



They are all equal and independent of n. 

In first order of the perturbation we have also non-zero off-diagonal matrix elements for transitions between different 
levels in the two groups. For n/mwe have^ 



(1)(12) = V 
k^n 



v (ll) ,(0)(12) 
nk J km 



v (22) ,(0)(12) 
km J nk 



V (H) f (0)(12) 
. Vnm J mm 



E, 



(0) 



E 



(0) 



k^m 



E. 



(o) 



E. 



(0) 



E 



(0) 



(0) 



(34) 



k—m 



k—n 



Here we have used the property 



y(U) = _y(22) = _ ^ (y \ 
Y nm nm D 1 V **» # 

hi \ / ? 



(35) 



which follows from the last term in Eq. I|32l) . 

The energy denominators in Eq. (|34|) correspond to the energy difference in one group of levels and are of the 
order of the small quadrupole energies AEq. In Eq. (|3*ljl we neglected the contributions of the off-diagonal matrix 

elements Vnm and the diagonal matrix elements /i!" 11 ', /n"^ 22 ^. This is justified since they have much larger energy 
denominators, equal to the distance between two levels from two different groups. This distance is of the order of 
the TLS energy E which is much larger than AEq. Therefore, this contribution is small due to the small parameter 
AE Q /E < 1. 

Taking into account that for real matrix I Vq ) = I Vq j , we get from Eq. I)34|l an important relation for the 

V / nm V / mn 

off-diagonal matrix elements with n m 



(1)(12) = 
nm 



-/, 



(1)(12) 
mn 



(36) 



(12) 
ij 



(12) r . , . 

-a) z for i f: j. 



It is equivalent to the important property that a 

In the second order in the perturbation we have also corrections to the off-diagonal transitions matrix elements 
/„' for transitions between equal levels in the two groups. Using second order perturbation theorjiiS we get after 
rather cumbersome calculations 



f(12) 
/ nn 



f(0)(12) 
/ nn 




(37) 



In this approximation the matrix elements / ; 
(El. 



(12) 



differ for different n and are smaller than the unperturbed values 



ECHO OSCILLATIONS IN A MAGNETIC FIELD 



In this section, using the above results, we will investigate the echo amplitude as a function of the time delay T21 
between the two electric pulses. This dependence straightforwardly give rise to the echo oscillations in the external 
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FIG. 4: A four-level system for nuclear spin J = 1/2. 



magnetic field. For ease of understanding, let us start from the simplest case of a four-level system. This corresponds 
to a nuclear spin J = 1/2 in a magnetic field. Such nuclei have zero quadrupole moment and, therefore, do not 
interact with a TLS. Nevertheless, formally we can consider this case to illustrate qualitatively the main features of 
our theory. We show that these features will be conserved in the more realistic case of J > 1/2. 



A. Four-Level System 

In the case of ./ = 1/2 we have the four-level system shown in Fig.0J Introducing the following notations for the 

( 12) 

off-diagonal matrix elements a ik between two group of levels (1) and (2) 



.,(12) 



,,(12) 



a (l 



2b 2 ) 



,,(12) 



,,(12) 



2a &, 



(38) 



where 



ao = e • m 



Ao 

E 1 



E e 



(39) 



and taking into account that in the framework of the perturbation theory b 1, we get from Eq. H15|) the following 
expression for the echo-amplitude from a four-level system 



P r 



echo 



c 



Aa 2 b 2 e i£T21 1 



Aa^b 2 coser 2 i |l 



Cat 



1 



oo!,2 • 4 £r 21 

326 sin 



0(b 4 ). 



(40) 



Here the coefficient C includes all those factors which are for the time being irrelevant. 

From Eq. I|40() we can see that the echo amplitude oscillates as a function of the time delay between the two pulses, 
T21. The frequency of these oscillations is determined by the interlevel splitting e between the levels 1 and 2. Since 
this splitting is a function of applied magnetic field the echo will experience oscillations when magnetic field changes. 

The average level of these oscillations is given by 



P ecbo = C4{l-l2b 2 ). 



(41) 



For sufficiently high magnetic fields, the interlevel splitting e oc H and therefore, according to Eq. I|39(l . b oc 1/H — > 
for H — > oo. This means that with increasing magnetic field the average echo amplitude increases and finally 
approaches the value of the TLS echo amplitude, Ccjq without a nuclear quadrupole moment. It looks as if a 
sufficiently high magnetic field switched off the TLS-quadrupole interaction. According to Eq. (|40|l we have the same 
value of the echo amplitude Ca^ in the limit ttl — > but in a final magnetic field. 



B. TLS-Quadrupole Multi-Level System 



In this section we will generalize our results to the case of an arbitrary TLS-quadrupole multi-level system. Using 
a notation similar to Eqs. 138J1 and 1|39JI 

a<W=ao | l-2£#„, ). : -m^/ = 'laj, (42) 
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where 



(43) 



" I LI I b 7— j J 

and e mn — E$) — En^ is the distance between the quadrupole levels (inside one level group!), we get for the echo 
amplitude 



P C c ho = CatU-^ ]T bl, 

\ n,m>n 



■ 4 S nm T21 \ .... 

sm — - — . (44) 



Here N = 2 J + 1 and the factor C, as before, includes all irrelevant details. The average (over T21) of the echo 
amplitude is given by 



echo 




(45) 



The quadrupole level splittings s mn are functions of the external magnetic field H. In the absence of any level 
degeneracy, in small magnetic fields when the Zeeman energy fxH <C AEq , these splittings acquire a small quadratic 
corrections of the order of (^H) 2 / AEq. Therefore, in this case the echo amplitude attains, as a function of the applied 
magnetic field, oscillations with a frequency of the order of /j, 2 Ht2i/?iAEq. In the case of a level degeneracy at zero 
magnetic field, the splittings of the originally degenerate levels are of the order /iH and, therefore, the frequency of 
the echo oscillations is of the order of HT21/T1. A simple estimate shows that the order of the magnitude of the period 
of these oscillations is in a full agreement with the experimental datai^. 

With a further increase of the magnetic field, the Zeeman energy, /j,H, becomes of the order of or bigger than the 
typical quadrupole energy splitting AEq in zero magnetic field. In such magnetic fields the typical level splittings 
£nm are mainly determined by the Zeeman energy fiH and in the limit fiH 3> AEq they become linear functions of 
H. Therefore, in this case the matrix elements b nm cx l/e nm ocl/H are dying-off with increasing magnetic field. As 
a result the average echo amplitude increases, approaching for H — > 00 the level of the TLS echo amplitude with zero 
quadrupole moment. 

Above we have considered the case of one nucleus interacting with a TLS. However it can happen in glasses that 
many nuclei participate in a tunneling motion. In such a case, since the magnetic interaction between the nuclei can 
be neglected, we can easily take this into account by including in Eq. H44fl a summation over the nuclei 

Pecho = C4(l-£^ £ (^jW^pV (46) 
\ s n.m>n / 

where N s = 2J S + 1. For distant nuclei the matrix elements b s nm depend on the distance, R s , between the nucleus 
and the TLS as b s nm oc 1/-Rg. Therefore, the main contributions to the sum over nuclei in Eq. 146JI are from nuclei 
participating in the tunneling and nearest neighbors. 



C. Difference Between Integer And Half-Integer Spins 

Up to now we did not discriminate between the two cases of integer and half-integer nuclear spin J. There is, 
however, an important difference between their quadrupole energy level spectra in zero magnetic field. In the case of 

an integer spin, J = 1, 2, the energy levels of the quadrupole Hamiltonian Wq (or Wq) in an electric field gradient 
tensor ipik(0) of arbitrary symmetry are not degenerate, whereas in the case of a half-integer spin, J = 3/2,5/2, 
according to Kramer's theorem all energy levels of the quadrupole Hamiltonian Wq are double degenerate. This 
degeneracy can only be lifted by applying a magnetic field. 

To calculate the echo amplitude we have used a perturbation approach for a non-degenerate case. Strictly speaking 
this approach is only valid in a finite magnetic field. Therefore, we now discuss the special case of a half-integer nuclear 
spin and H = 04i In this case the energy denominators of matrix elements fnm 12 ^ and b nm (Eqs. Q34J1 and (@3|) are 
zero if n and m belong to a pair of Kramer's degenerate levels. On the other hand the matrix elements I Vq ) are 

also zero in this case. The reason for this is obvious. If the matrix elements were not zero the perturbation Vq (or 



10 



Vq) would lift the Kramer's degeneracy of the nuclear spin levels given by Wq. But this cannot be the case since the 
interaction of the quadrupolc moment with a TLS motion Vq belongs to the same symmetry class as Wq. Therefore, 
the matrix elements belonging to a pair of Kramer's degenerated levels n and n' for H = one has for H = obey 
the relations 

s(o)\ /s(o)\ /s(o)\ 

V Q ) =0 for n ^ n', and [V Q =[V Q ) . (47) 



As a result, we have in our equation for half-integer spin and H — an indeterminate form 0/0 which should be 
properly evaluated. For this we need the limiting behavior of both, numerator and the denominator, as H — * 0. It is 
clear that in a general the denominator, being the energy splitting between two Kramer's levels in a magnetic field, 
vanishes proportional to H . The same is true also for the numerator. The analysis, based on a perturbation theory for 
degenerate levels^ gives the following expression for non-diagonal matrix elements describing the transition between 
two Kramer's degenerate levels, n, n' , split in an arbitrarily small magnetic field H 

^ (n H ) nm (v^) , + (^ 0) ) (n H ) mn , 

/ in 1 -) \ \ \ '_ mn' \ _ / nm (AR\ 

[ V Q ) nn , ~ 2. E W -E W ' ( } 

m^n,n> a n ^m 

Here the superscripts (0) and (1) refer to the matrix elements in zero and non-zero small magnetic fields H, respectively. 
(Hff) nm are the matrix elements of the Zeeman term, Eq. l(T8jl, 

We see from Eq. (|48fl that the matrix elements ( V^p ) are indeed proportional to the magnetic field H at small 

V w J nn' 

fields. Therefore, we come to the important conclusion that at zero magnetic field the ratio ( V Q ) /s nn i has a 

\ W J nn' 

finite value which depends on the orientation of the magnetic field H relative to the principal axes of the internal 
electric field gradient tensor f^ifc(O). Thus, for H = 0, the coefficients b nn > for the transitions between Kramer's 
degenerate levels are also finite and depend on the orientation of the applied magnetic field H. In other words, we 
find in this case a non-analytical behavior of the functions &„„/(H) at H — > 0. 

This behavior, for half-integer spin, of the matrix elements b nn i in a weak magnetic fields has important conse- 
quences. For small magnetic fields, when the Zeeman energy fiH is much smaller than the quadrupole energy AEq, 
the distance between the split originally degenerate levels is of the order of <§; AEq. On the other hand the 
distance between other (non-degenerate) levels is of the order of AEq . According to Eq. (|44|l a contribution to the 
echo signal from a transition n — > m is proportional to the oscillating function (— b 2 nm ) sin (s nm T2i/2). This function 
has a maximum for e nm — 0. Therefore, for small magnetic fields, the contributions of all quasi-degenerate pairs of 
levels are in phase and proportional to — sin 4 (/iiiT2i/2). With increasing magnetic field this function decreases from 
the value to the minimum value -1 for \iHti\ = 7r. In our view, this explains why in experiment there is often a 
sharp maximum in the echo amplitude at zero magnetic fieldi^. We believe that this maximum is a contribution 
from Kramer's degenerate pairs of levels split by weak magnetic fields and indicates the presence of nuclei with a 
half-integer spin J in the glass. 



VI. FINAL REMARKS 



In this final section we address shortly the main approximations we have made and the limitations of the presented 
theory. As already mentioned in the beginning of the paper we have skipped in our equations all resonance factors like 
ij^l or This is possible when all such factors are similar for all transitions between two group of levels. For this, 
the spectral width h/n of the electric pulses should be much larger than the typical quadrupole energy splitting AEq 
which means that the electric pulses should be sufficiently short, AEQTi/K <ti 1. Otherwise the important property 

(12) (12) 

<X-j — for i 7^ j which we used in the paper is not sufficient. The reason is that i — ► j and j i transitions 

(121 ('121 

correspond to different energy differences (see Fig. 0] for example) and one therefore should multiply a\j and 
with different resonance factors (3ij. In such a case the final expression for the echo amplitude becomes much more 
cumbersome. 

The other approximation we made was that we considered the limit of a rather small electric field amplitudes when 
the Rabi frequencies (F1.2 -m) Aq/TlE are much smaller than ft/r^. In the case of short electric pulses (as was indicated 
above) this limitation can be easily overcome and one can show that Eq. I|44|l for the echo amplitude will remain valid 
and only the coefficient C changes. 
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